Threshold of coexistence and critical behavior of a predator-prey cellular automaton 
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We study a probabilistic cellular automaton to describe two population biology problems: the 
threshold of species coexistence in a predator-prey system and the spreading of an epidemic in a 
population. By carrying out time-dependent simulations we obtain the dynamic critical exponents 
and the phase boundaries (thresholds) related to the transition between an active state, where prey 
and predators present a stable coexistence, and a prey absorbing state. The estimates for the critical 
exponents show that the transition belongs to the directed percolation universality class. In the limit 
where the cellular automaton maps into a model for the spreading of an epidemic with immunization 
we observe a crossover from directed percolation class to the dynamic percolation class. Patterns of 
growing clusters related to species coexistence and spreading of epidemic are shown and discussed. 

PACS numbers: 05.70.Ln, 87.23.Cc, 64.60.Ht 
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I. INTRODUCTION 

In 1958, Huffaker 1], in a pioneering experiment, has 
been able to maintain in the laboratory a population 
of prey and predators coexisting and presenting self- 
sustained coupled time oscillations. He has verified that 
persistence of species was only possible in a large and het- 
erogeneous space. Since then different models have been 
proposed to explain the role of space in determi ning the 
species coexistence jllSlllllSlllilillSinilll Ellli 
A common feature of these models is that they are based 
either on interacting particle systems 0, 0| , also called 
stochastic lattice models 0, or on probabilistic cellular 
automata. These descriptions are appropriate and rele- 
vant when considering predator-prey systems which are 
under conditions of very low species population densities 
and/or when their habitat is spatially heterogeneous 2]. 

In the present article we study a stochastic model 
that takes into account the spatial structure explicitly 
to describe two population biology issues: the thresh- 
old of species coexistence in a predator-prey system and 
the spreading of an infectious disease in a population. 
The individuals of each population are treated as dis- 
crete and they are supposed to occupy the sites of a 
two-dimensional lattice. The interactions between indi- 
viduals are included in the local stochastic rules which 
take into account the Lotka-Volterra mechanisms of in- 
teraction 0, 01 . The dynamics associated to the sys- 
tem is a Markovian discrete time process defined by a 
probabilistic cellular automaton, having three states per 
site, which we call predator-prey cellular automaton. The 
connection between pattern formation and species coex- 
istence in this model was recently analyzed [iJI ■ Results 
from steady- state simulations have shown that, depend- 
ing on the set of the parameters of the model, the follow- 
ing steady states can be attained 0|: a prey absorbing 
state, where predators have been extinct and the lattice 
is full of prey, and an active state where both species co- 
exist with constant time densities. For finite systems, it 



was also detected an active state where both population 
densities oscillate in time with the same characteristic 
frequency pTl |. 

Our aim in this work is to determine with precision the 
threshold of stable coexistence of species in the predator- 
prey cellular automaton. To this end we perform time- 
dependent simulations [HIMilllllilillillillil 
l28i . .291 . .30. .31.] which allow us to obtain the phase bound- 
aries together with the critical exponents. The scaling 
analysis of the time-dependent simulations yields a set 
of dynamic critical exponents and thus the possibility of 
classification of the models with absorbing states in uni- 
versality classes. Here, we have found that the automa- 
ton exhibits a line of continuous phase transition which 
belongs to the universality class of directed percolation, 
for all sets of parameters such that the prey birth prob- 
ability is different from zero. When this probability van- 
ishes there occurs a crossover to the universality class of 
dynamic isotropic percolation. In this limit the predator- 
prey cellular automaton is mapped into a model for the 
spreading of an epidemic with immunization, a general 
epidemic process |i20i, i21|, jS^ 33, 34]. 

This article is organized as follows. In Sec. II we 
describe the predator-prey cellular automaton and the 
time-dependent simulations. Critical properties and pic- 
tures of growing clusters generated by the simulations are 
shown in the same section. In Sec. HI we analyze the 
critical behavior of a model for an epidemic spreading. 
We briefly discuss and summarize our results in Sec. IV. 



II. PREDATOR-PREY PROBABILISTIC 
CELLULAR AUTOMATON 

A. The model 

We assume that each individual of each species popu- 
lation can reside on the sites of a regular square lattice 
which represents their habitat. A site in the lattice can 
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be in one of three states: occupied by a prey individ- 
ual (X); occupied by a predator individual (Y) or empty 
(Z). The predator-prey probabilistic cellular automaton 
comprehends the following processes: 

Z + X^2X, (1) 

X + Y ->2Y, (2) 

and 

Y ^ Z. (3) 

The transitions between states obey probabilistic rules 
which depend on the state of the given site and on the 
states of its four nearest neighbors at the north, east, 
south and west (which defines the neighborhood). The 
update is synchronous and the rules are: 

(i) A prey individual can be born in an empty site with 
a probability a/4 times the number of sites occupied by 
prey in the neighborhood, process 1^. 

(ii) A predator individual can be born in a site occu- 
pied by a prey if there are predators in its neighborhood. 
Prey disappear instantaneously and give place to a new 
predator. The probability of this process is 6/4 times the 
number of sites occupied by predators in the neighbor- 
hood, process 

(iii) The death of predators is spontaneous: a site oc- 
cupied by a predator can be evacuated with a probability 
c. This process reintegrates to the system the resources 
for prey prohferation, process ©. 

The model has three parameters a, b, and c, with 
0<a<l, 0<&<landO<c<l. However, we 
assume that a + b + c = 1, thus just two parameters are 
independent. We consider the parametrization 0, [l3 |: 
a — (1 — c)/2 — p and b = {1 — c)/2 + p. The param- 
eter p is such that —1/2 < p < 1/2 and parameter c 
denotes the death probabihty of predators, < c < 1. 
This parametrization allows us to analyze the model in 
a triangular phase diagram p — c, as shown in Fig. ^ 

The system evolves in time and can eventually reach 
stationary states characterized by constant time densities 
of prey, predators, and empty sites. These can be either 
active states where prey and predators have stable coex- 
istence or absorbing states. We note that finite systems 
also exhibit local oscillatory behavior of the densities of 
prey and predators [Tll |. 

A transition line ci (p) crosses the entire p — c diagram 
represented in Fig. ^ It starts from the left corner of the 
triangle and ends on the opposite side. The densities of 
prey, predators and empty sites change continuously at 
ci{p). For a ^ 0, the transition is from the active state 
to the prey absorbing state and the transition line was 
estimated firstly in reference 0|. In the limiting case 
corresponding to a = 0, the point where ci{p) meets the 
right side of triangle of Fig. the transition is from 
one of the infinitely many absorbing states to the prey 
absorbing state. In order to determine the critical line 
ci(p), with precision, and the dynamic critical exponents, 
we perform time-dependent simulations. 
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FIG. 1: Phase diagram in the plane p — c, showing regions cor- 
responding to the prey absorbing phase and the active phase. 
Phases are separated by a transition line obtained from time- 
dependent simulations (see Table I). 



B. Critical exponents and thresholds 

To carry out a time-dependent simulation analysis for 
the predator-prey cellular automaton we follow the time 
evolution of states very close to the prey absorbing state. 
We depart from an initial condition (t = 0) with a single 
predator at the origin of a lattice covered by prey. Once 
the system is placed in the initial condition we apply 
the local rules (i), (ii), and (iii), described in Sec. II A, 
and let it evolves in time. We have considered Ng — 
10000 samples (independent runs) all starting with this 
same condition. For each fixed value of p we vary c near 
its transition value. The simulations were performed in 
systems sufficiently large so that predators do not reach 
the borders of the lattice. 

We have investigated the behavior of the following 
quantities: the survival probability of predators P, that 
is, the probability that predators have not been extinct 
until time t; the mean-number of predators (A^) at time 
t, whose average is calculated over all the samples, in- 
cluding those where predators have been extinct before 
time t\ and the average mean-square distance of spread- 
ing of predators from the origin, (i?^). This average is 
calculated by taking into account only the samples that 
survived until time t. 

According to the scaling laws for time-dependent simu- 
lations 1!^, for large values of t, we expect, at the critical 
point, the following power laws, 

(TV) r^t\ (4) 

P ^ t-\ (5) 
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and 



(6) 



where rj, S and z are the critical exponents associated to 
the mean-number of predators, the survival probability, 
and the mean-square distance of spreading, respectively. 

As can be seen in Fig. El the log-log plot of the mean- 
number of predators, as well as of the survival probability, 
show very well defined critical and off-critical asymptotic 
behaviors. For a fixed value of p, and values of c below, 
but close to the critical value (supercritical regime), these 
quantities present a positive curvature, which indicates 
the presence of activity in the system; and for values 
of c greater than the critical value (subcritical regime) 
they present a negative curvature. At criticality, (N) and 
P exhibit a very clear asymptotic power law behavior. 
Therefore, it is possible to obtain, with precision, the 
phase boundaries together with the critical exponents. 
From Fig. |21we get the following estimates for the critical 
exponents 



r] = 0.230(2) 



(5 = 0.451(3), z = 1.134(4), 



for p = (that is, for a ~ b). As reported in Table I, 
the same values for the exponents, with the error bars, 
were obtained (using the same procedure) for different 
values of p along the critical line ci{p) (excluding the 
critical point for which a = 0). The errors of the ex- 
ponents were estimated by dividing Ns — 10000 sam- 
ples in bins of 2000 runs, and calculating the fluctuation 
of the averages obtained for each bin. These values of 
the dynamic critical exponents are consistent with the 
ones of the directed percolation class in 2 -f 1 dimensions: 
r] = 0.2295(10), S = 0.4505(10) and z = 1.1325(10) H. 

This dynamic critical behavior has been obtained for a 
diversity of models with finite number of asymmetric ab- 
sorbing states ,27>.29,j . such as the contact process, which 
can be viewed as a model for a simple epidemic |15). 
The critical behavior of this model has been extensively 
studied and it belongs to the directed percolation class 
|26l|27|,l28|. Other examples are the Domany-Kinzel cel- 
lular automaton |35l and the spreading of damage tran- 
sitions [27113^ l37|. Among the models with more then 
one absorbing state, per hap s the most known is the Ziff- 
GuUari-Barshad model jS^I for the reaction of oxidation 
of carbon monoxide over a catalytic surface; it also has 
dynamic critical exponents l^lj consistent with the ones 
associated to the directed percolation class. 

The critical exponent associated to the time correla- 
tion length can be obtained from the time-dependent 
simulations. Here, the estimation was accomplished by 
studying the time behavior of the derivative D of the 
mean- value of predators (N). From the scaling laws it 
can be shown that Hal , 



d\og(N) 

D = — — - ~ i / II 

dlogc 



(7) 
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at the critical value of c. Data for different values of c, 
close to the critical point, were taken in the same run 



FIG. 2: Log-log plot of the mean-number of predators (A'') 
(upper panel), the survival probability P (middle panel) and 
the mean- square distance of spreading of predators (_R^) 
(lower panel), as a function of the time t, for p = 0. Each 
figure shows the behavior of these quantities for diflterent val- 
ues of c. From top to bottom: c = 0.18875, c = 0.18975, 
c = 0.19075, c = 0.19175 and c = 0.19275. 




FIG. 3: Logarithmic derivative D of (TV) with respect to c. 
For p = and considering values of c slightly above and below 
the critical value. 



TABLE L Phase boundaries and dynamic critical exponents 
for the predator-prey cellular automaton. The values of the 
last column were obtained using 13 = Suu . 
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1. 



13 



-0.495 

-0.45 

-0.40 

-0.30 

-0.20 

-0.10 

0.00 

0.10 

0.20 

0.22 

0.24 

0.26 

0.28 

0.30 

0.32 

0.34 

0.36 

0.38 



0.00206(3) 

0.0228(3) 

0.04515(10) 

0.0877(2) 

0.12695(2) 

0.1617(2) 

0.19075(10) 

0.2127(2) 

0.2243(2) 

0.2253(2) 

0.2262(2) 

0.2263(2) 

0.2262(2) 

0.2257(1) 

0.2250(2) 

0.2239(2) 

0.223(1) 

0.223(2) 



0.234(4; 
0.227(4; 
0.229(3 
0.226(3 
0.232(2 
0.234(3 
0.230(2 
0.230(2 
0.228(2 
0.231(2 
0.230(3 
0.226(4 
0.232(4 
0.231(4; 
0.233(4; 
0.234(5 
0.235(7; 
0.24(1) 



0.444(6 
0.448(4 
0.451(3 
0.452(3 
0.452(3 
0.447(3 
0.451(3 
0.449(3 
0.450(3 
0.447(3 
0.445(5 
0.443(6 
0.444(6 
0.445(5 
0.448(5 
0.449(6 
0.452(8 
0.44(1) 



1.13(2) 

1.13(1) 

1.135(8 

1.134(7 

1.129(8 

1.135(6 

1.134(4 

1.131(7; 

1.128(8 

1.131(7; 

1.136(9 

1.129(8 

1.135(10) 

1.13(1) 

1.14(1) 

1.13(1) 

1.14(2) 

1.14(2) 



1.30(2) 

1.31(2) 

1.30(1) 

1.29(1) 

1.288(9) 

1.295(10) 

1.292(8) 

1.287(8) 

1.291(9) 

1.29(1) 

1.30(1) 

1.29(1) 

1.29(1) 

1.295(10) 

1.30(1) 

1.29(1) 

1.28(2) 

1.29(2) 



0.58(2) 

0.59(1) 

0.586(8) 

0.583(8) 

0.582(8) 

0.579(8) 

0.583(7) 

0.578(7) 

0.581(8) 

0.58(1) 

0.58(1) 

0.57(1) 

0.57(1) 

0.58(1) 

0.58(1) 

0.58(1) 

0.58(2) 

0.57(2) 



0.00 0.220(3) 0.587(2) 0.096(4) 1.767(5) 1.51(1) 



0.145(7) 
0.139(4)t 



tSee Sec. Ill D. 



and the derivatives calculated considering finite differ- 
ences; the behavior of for p = 0, is shown in Fig. 13 
As reported in Table I, the value obtained for the expo- 
nent along the transition line, is in agreement with 
the estimated value ^|| = 1.295(6) jl^ for the directed 
percolation in 2 -I- 1 dimensions. 

The static exponent (3 associated to the order parame- 



FIG. 4: Two snapshots of the predator-prey automaton, gen- 
erated from a single predator at the origin (center) of a lattice 
covered by prey, for p — —0.45 and c — 0.01 (upper panel, 
taken at t = 10000) and c = 0.0227 (lower panel, taken at 
t = 65000). The white points represent sites occupied by 
prey, the red points by predators and the blue points are 
empty sites. 



ter was calculated from the values for the dynamic critical 
exponents S and and using the scaling relation P = Sv\\ 
19, 27|]. The resulting value for the exponent f3, reported 
on Table I, is consistent with directed percolation in 2-1- 1 
dimensions p^ . 

For sets of parameters such that a approaches zero 
the critical behavior of the model suggests a crossover to 
another universality class. This feature will be explored 
in Sec. III. 



C. Growing clusters of predator-prey coexistence 

The pictures shown in Figs. andlHlare related to 
configurations generated by simulations of the predator- 
prey cellular automaton departing from the initial con- 
dition with one predator at the origin and the lattice full 
of prey. It can be seen in Fig. ^ a configuration in the 
supercritical regime (upper panel) and an almost critical 
configuration (lower panel) for a set of parameters which 
corresponds to 6 << a and c << a (left down corner of 
the triangular phase diagram of Fig. In this case. 




FIG. 5: The same of Fig. 4, for p = -0.3 and c = 0.0075 
(upper panel, taken at t = 1750), and c = 0.0867 (lower 
panel, taken after t = 20000). 



predators stay in the lattice for a long period of time, 
and when one of them dies it gives place to an empty site 
that will be almost immediately occupied by a prey indi- 
vidual. Consequently, just a few number of empty sites 
are present in the steady state and also in the growing 
clusters. This behavior suggests that the present three 
state process, in the limit b —>■ 0^ c —> with b/c finite, 
can be replaced by a two state process, similar to the 
contact process, involving predators and prey, as pointed 
out in reference [l^ . 

The growing cluster of the upper panel of Fig. [S] shows 
a pattern formation related to an active region where 
prey and predators coexist and can exhibit local time 
coupled oscillations. With a nonzero probability, the 
dynamics will survive forever with the densities of each 
species different from zero but less then 1. This cluster, 
as well as the supercritical growing clusters of Figs, ^and 
EK upper panels), assumes after some time, which we call 
T, the asymptotic shape of a ball. The correlation length 
^ can be understood as the linear size of the cluster at 
t = T so that for t > t the system is close to the station- 
ary active state. Our results suggest that the radius R 
of the ball grows linearly with t for t > r in accordance 
with the shape theorem j^. At the critical point, t and 
^ diverges. Far from the critical point, the time t the 
cluster takes to get a defined shape can be very small. 
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FIG. 6: The same of Fig. 4, for p = 0.3 and c = 0.06 (upper 
panel, taken at t = 500), and c = 0.2257 (lower panel, taken 
at t = 3000). 



For instance, for the set of parameters considered in the 
upper panel of Fig. El the shape of a ball is already taken 
after just 200 time steps. 

We notice that the supercritical cluster of Fig. cor- 
responds to a configuration where the populations of 
prey and predators, under conditions of low densities, 
are grouped into small clusters of a unique species that 
are isolated from each other, and coexist without extinc- 
tion in a highly heterogeneous space. In this respect, this 
species coexistence can be compared to the Huffaker's ex- 
periment, commented in Sec. I. The coexistence for the 
set of parameters of Fig. El is associated to local self- 
sustained coupled time oscillations of prey and predators 
populations and has been described in detail in reference 

Typical clusters at criticality are shown in the lower 
panels of Figs. 01 and |S1 According to our results 
of Sec. II A, these configurations must be related with 
transitions belonging to the universality class of directed 
percolation. They do not present a ball shape, as the su- 
percritical clusters described above, but rather a shape 
of a fractal nature. We note that the critical cluster of 
Fig. El (lower panel) presents, in contrast with those of 
Figs. ^ and El (lower panels), noticeable agglomerations 
of empty sites. For values of p corresponding to small val- 
ues of a, and, therefore, inside the crossover region, the 
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FIG. 7: Two snapshots of the epidemic with immunization 
(a = 0) generated from a single infected individual located 
at the origin (center) of a lattice covered by susceptible indi- 
viduals, for c = 0.15 (upper panel, taken ai t — 500) and for 
c = 0.22 (lower panel, taken ai t = 1000). The white, red and 
blue points represent sites occupied by susceptible, infected 
and an immune individuals, respectively. 



critical clusters present larger agglomerations of empty 
sites. In the first time steps of the simulation, they re- 
semble the critical growing cluster for a = (lower panel 
of Fig. [TJ. However, as long as the parameter a is differ- 
ent from zero, even if it is very small, the critical cluster 
will present agglomerations of empty sites which can be 
always active, inasmuch as prey can always proliferate 
in empty sites. We remark that these critical clusters 
are still related to transitions belonging to the directed 
percolation class. 



A qualitative similar behavior of critical clusters was 
obtained by Dammer and Hinrichsen |34| for an epidemic 
spreading model which also presents a crossover from the 
directed percolation critical behavior to the dynamic per- 
colation critical behavior. 



III. SPREADING OF AN EPIDEMIC 
A. The model 

For a = the predator-prey probabilistic cellular au- 
tomaton can be interpreted as a model for the propaga- 
tion of an epidemic with immunization. It mimics the 
spreading of an epidemic in a population composed by 
susceptible individuals {X) that become infected (Y) by 
contact with infected individuals; once infected the in- 
dividuals can recover, and become immune (Z) sponta- 
neously. The process of infection can be represented by 
the reaction: 

X + Y ^2Y, (8) 

and the recovery process, which includes immvmization, 
by the reaction: 

Y^Z. (9) 

These are the basic processes which are taken into ac- 
count in the modeling of the spreading of an epidemics 
with immunization [53, |32, 113; 113 • The present 
model is defined on a regular square lattice where each 
site can be in the following states: occupied by a suscepti- 
ble, an infected or an immune individual. It comprehends 
the following stochastic rules: 

(A) The infection can occurs when a susceptible indi- 
vidual, which occupy a given site, has at least one site 
occupied by an infected individual in its neighborhood, 
reaction This process occurs with probability 6/4 
times the number of infected individuals in the neighbor- 
ing sites. 

(B) The recovering process can occurs spontaneously 
with probability c when a site is occupied by an infected 
individual, reaction ©. The condition b + c = 1 is 
obeyed, with b being the infection probability and c the 
recovery probability. This model can exhibits infinitely 
many absorbing states and presents a continuous tran- 
sition belonging to the dynamic percolation universality 
class, as we show in Sec. Ill C. 

B. Patterns of spreading 

To characterize the threshold of the epidemic spread- 
ing we perform time-dependent simulations. Initially all 
sites of the lattice are occupied by susceptible individuals 
except one, in the origin (center) of the lattice, which is 
occupied by an infected individual. The stochastic dy- 
namics follows the rules (A) and (B) defined in Sec. Ill 
A with a synchronous updating. The infection starts to 
be transmitted at time t — 0. The problem is to know 
how an epidemic propagates (or not) in the population. 
For high values of the probability b of infection compared 
to the values of the recovering probability c, the epidemic 
spreads leaving a cluster of inactive sites composed by im- 
mune individuals and some groups of individuals which 
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FIG. 8: Log-log plot of the mean-number of infected (up- 
per panel), the survival probability (middle panel) and the 
mean-square distance of spreading of infected (lower panel), 
for a = 0. Each figure shows the behavior of the quantities for 
different values of c. From top to bottom: c = 0.210, 0.215, 
0.220, 0.225 and 0.230. 
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FIG. 9: Logarithmic derivative D of (TV) with respect to c. 
For a = and considering values of c close to the critical 
value. 

remain susceptible forever. A typical picture of the lat- 
tice in this situation in shown in the upper panel of Fig. 
□ The cluster grows with a front of infected individu- 
als which stay in the border. After a finite time r, the 
cluster assumes a limiting shape of a ball which spreads 
to infinity with a nonzero probability. We observe that 
different initial seeds lead to different configurations of 
the spreading of the epidemic, giving rise to an infinitely 
many absorbing states. As h is decreased (and c is in- 
creased) the threshold for the spreading of the epidemics 
is reached. Above the threshold, the epidemic will cease 
in a finite time leaving a cluster with just a few number of 
immune individuals, the rest of the lattice being covered 
by susceptible individuals. A picture of the lattice near 
the threshold of epidemic spreading is shown in the lower 
panel of Fig. [7| This almost critical cluster presents an 
irregular shape of fractal nature. 

Similar critical and supercritical clusters have been ob- 
tained from time-dependent simulations for other spatial- 
structured stochastic models for an epidemic with immu- 
nization [lll27Llp|. 

C. Dynamic critical exponents 

In Fig. |S1 it is shown the time behavior of the mean- 
number of infected individuals {N) , the survival probabil- 
ity P, and the mean-square distance of spreading from 
the origin {B?) obtained from time-dependent simula- 
tions performed near the critical point. Since the log- log 
plots of of these quantities show very clear power-law be- 
havior, associated to critical behavior, then the dynamic 
critical exponents 77, 5 and z, defined in Eqs. Q), © and 
can be obtained with precision. The estimated crit- 
ical exponents are reported in Table I. These exponents 
are consistent with the critical exponents of the dynamic 
percolation universality class [29|. This is the expected 
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FIG. 10: Density of immune individuals po as a function of c, 
for a = 0.0. From top to bottom (at right): L = 10, 20, 30, 40, 
60, 80, 120, 160, 240. The critical value of c for the infinite 
system is Cc = 0.220(3) (see Table 1). 
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FIG. 11: Log-log plot of po, calculated at the critical value 
Cc = 0.220, for a = 0.0, as a function of L. Simulations were 
performed on square lattices of size ranging from L = 10 until 
320, and considering periodic boundary conditions. 



percolation [23|. We also have estimated this exponent 
by performing steady state simulations in finite systems, 
as described in the next section. 



D. Critical exponent /9 

To obtain another estimation of the critical exponent 
P associated to the order parameter we have considered 
square lattices with linear size L and periodic boundary 
conditions. For each value of L we have performed several 
independent simulations runs all starting from a lattice 
covered with susceptible individuals with the exception 
of one infected at the center of the lattice. The system 
evolves in time according to the synchronous stochastic 
dynamics defined in Sec. Ill A. In principle, any one of 
the infinitely many absorbing states can be reached; and, 
consequently, the number of immune individuals at the 
steady state varies from sample to sample. A simulation 
is finished when the system enters in an absorbing state 
and then the number of immune individuals is calculated. 
The mean value of this quantity, divided by the total 
number of lattice sites L^, is the density po of immune 
individuals at the steady state (for a given lattice size). 

The behavior of po as a function of c, for different val- 
ues of L is shown in Fig. ^| In the subcritical regime (c 
above its critical value Cc) po decreases as L is increased 
and in the limit L — > oo it vanishes. In the supercritical 
regime po is almost independent of L. We have assumed 
that Po is an appropriate order parameter for this transi- 
tion. Assuming the finite size scaling hypothesis (2^ we 
expect that po(A,L), where A = c — Cc, calculated for 
each value of L, scales at the critical point as, 

PoiO,L)^L-^/''\ (10) 

where is the critical exponent related to the spatial 
correlation length. 

To obtain the ratio /3/i^± we plot logpo, calculated 
at Cc, versus \ogL, as shown in Fig. 1111 We obtain 
l3/i^i_ = 0.185(5). Using the result = 4/3 ^ we find 
the value /3 = 0.139(4) which is in agreement with the 
one for the isotropic percolation, namely /? = 5/36 [2^ . 
And, also in agreement with the value estimated in Sec. 
Ill C. 



critical universal behavior for a general epidemic with 
immunization j20ll2l| . 

The critical exponent was obtained by analyzing 
the behavior of the derivative D, defined in Eq. ||7J), 
at the critical point. From Fig. |51 we obtain the esti- 
mation I'll = 1.51(1) which is consistent with the value 
I'll = 1.506(1) 29] of the dynamic percolation in two di- 
mensions. The static critical exponent f3 associated to 
the order parameter can be evaluated from the estimates 
for the dynamic critical exponents using the scaling re- 
lation P — (Ji'ii . The resulting value, which is given in 
Table I, is in agreement with the one for the dynamic 



IV. SUMMARY 

We have studied a spatial-structured model in which 
prey and predators individuals reside on the sites of a 
lattice and are described by discrete dynamic variables. 
The system evolves in time according to a probabilistic 
cellular automaton which takes into account the Lotka- 
Volterra interactions by the use of Markovian local rules. 
The model exhibits an active phase where prey and 
predator coexist without extinction, an absorbing prey 
phase, and, when the birth probability of prey vanishes, a 
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phase with infinitely many absorbing states composed by 
empty (immune) sites and prey (susceptible) . In the last 
case, it is more appropriate to refer to the predator-prey 
cellular automaton as a model for the spreading of an 
epidemic with immunization in a population composed 
by susceptible, infected and immune individuals. The 
automaton was studied by time-dependent simulations 
which allow us to determine with precision the phase 
boundaries (thresholds) and the dynamic critical expo- 
nents. We have localized a transition line which crosses 
the entire phase diagram. For all sets of parameters such 
that the prey birth probability is different from zero, this 
line separates the active phase, where prey and predators 
coexist, and the prey absorbing phase. We have shown 
that this transition belongs to the directed percolation 
universality class. We also have shown that when the 
prey birth probability equals zero, there is a crossover to 



the universality class of the dynamic percolation. 

Patterns of growing clusters with coexistence of prey, 
predators and empty sites (or, susceptible, infected 
and immune individuals), generated by the time depen- 
dent simulations, in the critical and in the supercritical 
regimes, were shown and discussed. The present study 
provides a detailed description of the thresholds of stable 
coexistence of the two species (or, the spreading of an 
epidemics) in the context of this predator-prey cellular 
automaton. 
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